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Abstract 

We analyze the decays of charginos, neutralinos, gluinos and the first/second gen- 
eration squarks in the Minimal Supersymmetric extension of the Standard Model, 
focusing on the three-body decays in scenarios where the ratio tan (3 of vacuum 
expectation values of the two Higgs doublet fields is large. We show that the three- 
body decays of the next-to-lightest neutralinos (lightest charginos) into hh, t^t^ 
{tv) final states, where third generation sfermion and Higgs boson exchange dia- 
grams play an important role, are dominant. Furthermore, we show that decays of 
gluinos into bb final states and squark decays into lighter sbottoms through gluino 
exchange can also have sizeable branching fractions, especially in scenarios where 
the soft SUSY breaking gaugino mass parameters are not unified at the GUT scale. 
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1. Introduction 



Supersymmetric theories (SUSY) |l| are the best motivated extensions of the Standard 
Model (SM) of the strong and electroweak interactions and the search for supersymmetric 
particles is one of the major goals of present and future collider experiments. The Minimal 
Supersymmetric extension of the Standard Model (MSSM), with R-parity conservation 
[leading to a stable lightest SUSY particle (LSP) which is in general the lightest neutralino 
Xi], the minimal particle content [three generations of fermions, two Higgs doublet fields 
to break the electroweak symmetry, as well as their SUSY partners] and a set of soft terms 
to break SUSY [e.g. mass terms for gauginos, sfermions and Higgs bosons, and trilinear 
couplings between squarks and Higgs bosons], is the most studied in the literature |T]]. 

In the MSSM, the phenomenology of the third generation sfermions is rather special. 
Indeed, due to the large value of the top, bottom and tau lepton masses, the two current 
eigenstates and Jr (with f = i,b or f) could strongly mix leading to a large 
splitting between the two mass eigenstates /i and /2 with the lighter one, /i, possibly 
much lighter than the other sfermions. In particular, for large values of the parameter 
tan P, the ratio of the vacuum expectation values of the two Higgs fieldsQ, and non-zero 
values of the Higgs-higgsino mass parameter fi, the mixing can be extremely strong in the 
sbottom and stau sectors, leading to relatively light 61 and fi states^. [In turn, for large 
values of the trilinear coupling At and/or small tan/5 values, the mixing can be strong 
in the stop sector, leading to a ti eigenstate much lighter than the other squarks and 
possibly lighter than the top quark itself.] 

A light bi [and/or fi] eigenstate may lead to dramatic consequences on the decay 
modes of the next-to-lightest neutralino, the lightest chargino, the gluino as well as the 
first and second generation squarks. Indeed, in this case, one has the following scenarii: 

{i) The possibility that squarks of the first and second generations decay into the 
lighter bi squark through the mode 

Ql,r ^ q g* ^ q bbl + qbbi (1) 

opens up. Because it is a QCD mediated decay, and because the decays into light charginos 
and neutralinos could be suppressed [e.g., if these particles are higgsino-like since the 
couplings are proportional to the light quark mass in this case], this mode might have 
a sizeable branching fraction, in particular in models where the gaugino masses are not 
unified at the GUT scale. 

(ii) If gluinos are lighter than squarks, they will mainly decay through virtual squark 
exchange into quarks and charginos/neutralinos; if bi is the lightest squark, its virtuality 
will be smaller, leading to the possible dominance of the decay mode 

g bb* ,bb bbx] and/or b t xf (2) 

^Note that large values of the parameter tan/3, tan/^ ~ 50, are favored in models with unification of 
the Yukawa couplings; see e.g. Ref. [Q. [The other solution, with tan/? ~ 1.5, seems to be excluded from 
the negative searches of MSSM Higgs bosons at LEP2.]. 

^Lighter third generation sfermions have been also advocated in many models such as those with 
inverted mass hierarchy where the first and second generation squarks are much heavier than their third 
generation partners; for recent papers see e.g. Ref. 



2 



(iii) The lightest chargino Xi aiid the next-to-hghtest neutrahno X2 "^ill decay into 
the LSP and two hght fermions. This occurs through gauge boson, Higgs boson and 
sfermion exchange diagrams. Because for high tan j3 values, the third generation sfermions, 
sbottoms and staus, are lighter and the Higgs bosons couple strongly to bottom quarks 
and tau leptons, the three-body decays: 

X'2 - Xibb, xU^r- (3) 
Xt - X?r+z/ (4) 

are in general enhanced compared to decays where first and second generation fermions 
are involved in the final states. For large values of tan/?, these three-body decay modes 
have been discussed in Ref. |^ and partly in Ref. 0] in the case of the neutralino X2- 

In this note, we analyze these various decay modes and investigate their consequences. 
In the next section we discuss the three body decays of squarks and gluinos, after pre- 
senting the analytical expressions for the partial widths. In Section 3, we analyze the 
three-body decays of the lightest chargino Xi the second hghtest neutralino x^'i 
give some illustrations on the branching ratios of these decay modes and compare our 
results with those of Ref. 0. Some conclusions are then given in section 4. 



2. Three— body decays of squarks and gluinos 

The Feynman diagrams for the three-body decays of squarks of the first and second 
generations into the lighter bi squark, eq.(l), and of gluinos through sbottom exchange, 
eq.(2), are given in Fig. 1. 





Figure: Generic Feynman diagrams for the three-body decays of squarks and gluinos. 

The Dalitz plot density of the decay mode eq. (1) is given in terms of the reduced 
energies of the two final state quarks, xi = 2Eq/mq. and X2 = 2Eb/mq-, and the gluino 
and sbottom reduced squared masses, fig = m?/m? and fi^ = m^/m'^_. Neglecting the 
masses of the final state quarksQ, and adding incoherently the contribution from the two 
different bbi and bb^ final states, it is given by: 

>.2 1 r _ 

Xi -X2 



dr(g, 



da;ida;2 



— rria. 

TT 



(1 - Xl - /ig)^ 



-xx-xi- - b\iig{\ -X1-X2- fil) 



(5) 



■^The expressions with the mass effects in this case and for gluino decays to be discussed later [with 
the additional final states g ^^xf etc.] are slightly more involved and will be given elsewhere 
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where the factors and 6j for the squark are given, in terms of the b mixing angle 6b, 
by: 



Oi = 62 = sin 6*6 and 02 = fei = cos 6^ 
Integrating over the energies xi, X2 with boundary conditions 

< xi < 1 — /i^ and 1 — — /x^ < X2 < 1 — ^J'l/i^ — xi) 
one obtains the partial decay width: 



(6) 



(7) 



3tt 



{fig - 3/i? + /ib + /i3/^b)Log 



l^g - 1 



2l4 



4fig 

{1 - flg){fll- fig) 



4 



+ ^{flg- fil + figfil)LOgfIl + ^ -{fig - 2/i? - 2fll + figfil) 



2flg 



Neglecting again the mass of the final bottom quark [but not in the couplings], the 
Dalitz plot density for the gluino decay g bbxi through 6j = &i,62 exchange is simply 
given by [here, Xi = 2Eh/mg,X2 = 2Ei/mg and /i^- = ml /m^, fx^ = m^Q/m^ 



dm 



xsuc ^ , 2 I t2 \ ^1 (1 ~ /^x ~ ^1) 



da;idx2 47r " V ' '''' ~ ^1 ^ f^k 

and after integration on the variables Xi and X2, one obtains for the partial width: 



(9) 



47r ^ 



+ (/^fe, - l)(^x - f^bjif^x - 3/^1 + f^k + /^b,/^x)Log^-^^ - /ijLog/i;^| (10) 
where the couplings between neutralino x?, bottom and sbottoms &i,&2, are given by 
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ce^ j V2Miycos/? 



;^j3 - 



V2 

— svi/(^ji — tan6'vK^j2) 



i3 



with = 1 — = sin^6'vi/ and sg^ = sin 6';,, etc.; the (rotated) matrix elements of the 
neutralino mass matrix are given in Ref. Note that for massless b quarks, there is 
no interference between the amplitudes with b\ and 62 exchange. 
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For illustration we analyze the branching ratios in two models, discussed in Ref. [10 
where the gaugino masses are not unified at the GUT scale: i) Models in which SUSY 
breaking occurs via an F-term that is not SU(5) singlet but belongs to a representation 
which appears in the symmetric product of two adjoints: (24(g)24)sym=l©24©75©200 
(where only 1 leads to universal masses); ii) The Oil model which is superstring motivated 
and where the SUSY breaking is moduli-dominated. The relative gaugino masses at the 
GUT scale and at the scale Mz are given in Tab. 1 taken from Ref. [PH]]. 





Ms 


Ms 


Ml 


1 


1(~ 6) 


1(~ 2) 


1(~1) 


24 


2(~ 12) 


-3( 6) 


-1(~-1) 


75 


1(~ 6) 


3(~ 6) 


-5( 5) 


200 


1(~ 6) 


2(~ 4) 


10(~ 10) 


on 


1(~ 6) 


5(~ 10) 


53/5(~ 53/5) 



Table 1: Relative gaugino masses at MqutCM^) in the F$ representations and the Oil model. 

In Fig. 2, we show the branching ratios of the decays qi qhh\ + qbhi (2a) and 
g -x^bb + x'^tb (2b) as a function of /i and for tan/3 = 50. In Fig. 2a, the (right- 
handed) squark mass is taken to be m^^ = 500 GeV while the gluino mass is slightly 
larger, rrig ~ M3 = 550 GeV. For large values of n, the lightest neutralinos and chargino 
are gaugino like, and in models 75 and Oil have masses comparable to the mass of the 
gluino [Mi,M2,M3 are of the same order in this case; the figure is cut in model Oil at 
fj, ~ M3, since the gluino becomes lighter than Xi]'- the two-body squark decays into 
charginos and neutralinos are phase space-suppressed and the virtuality of the gluino is 
not very large, leading to a rather large BR(g qbbi), reaching unity for the extreme 
values of /i where bi is light. In models 1 and 24, gluinos are much heavier than the lightest 
charginos and neutralinos [M3 3> Mi], and the two-body decays q ^ qx dominate [but 
BR(gi qbb\) reaches the level of a few percent]. Model 200 is an intermediate case. 
For small values of /x, the lightest chargino and neutralinos are higgsino like and much 
lighter than the gluino, 61 is also heavier and the two-body decay q ^ qx dominate. 

In Fig. 2b, the branching ratio g —>■ bbx^ + btx"^ is shown for M3 ~ rrig = 350 GeV 
and a common squark mass rriq = 600 GeV; with increasing /i the lightest sbottom mass 
varies from ~ 560 GeV (for = 200 GeV) to ~ 357 GeV (for /i = 1 TeV). Even for 
the universal gaugino mass case, the branching ratio is very large, especially for small 
and large values of /x. For /x ~ 200 GeV, 61 is not much lighter that the other squarks, 
but the x^bb coupling is enhanced [oc mb/Mvi/ tan /3; see eq. (11)] and the bb final state is 
favored. For large values of /i, m^^ becomes significantly smaller than rrig and the sbottom 
exchange channel dominates leading to the dominance of x^bb and x^tb final states. As 
in the case of squark decays, the situation is more favorable for bottom quark decays of 
the gluino in models 75 and Oil [in this case also we have cut the plot at /i ~ M3 since 
then, rrig < m^o], while it is similar for the 24 and 200 models. For larger fi values, /i > 1 

TeV, the two-body decay g —>■ bbi is accessible kinematically and will of course dominate 
[since it is a QCD process] all other channels. 
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Figure 2: The branching ratios BK{qR qhh\ + qhh) for a common scalar squark mass 
rriq = 500 GeV and M3 = 550 GeV (a) and BR(^ ^ x^fe + x^tb) for = 600 GeV 
and M3 = 350 GeV (b), as a function of /j, for tan /3 = 50 in the various models discussed 
above. M3 is fixed and the values of Mi and M2 are given Table 1. 
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3. Chargino and Neutralino Decays 



We have also calculated the partial decay widths of the three-body decays of the lightest 
chargino and of the next-to-lightest neutralino into the LSP and two fermions, eqs. (3-4). 
These decays are mediated by the exchange of gauge bosons [W boson for xt and Z boson 
for X2]) Higgs bosons [the charged boson for xt and the three neutral Higgs bosons 
h, H and A for as well as t and u channel sfermion exchanges [in particular the lighter 
fi slepton for xt and bi, fi states for Fig- 3. 



Figure 3: Generic Feynman diagrams contributing to the three-body decays of charginos 
and neutralinos into the LSP and two fermions. 

We have taken into account the contributions of all the channels [and of course the 
interference terms], the full dependence on the masses of the final fermions [to be able 
to describe more accurately the cases of chargino decays into rz/, tb and the neutralino 
decays into bb, ti final states] and the mixing of the third generation sfermions. Complete 
analytical formulae for the Dalitz plot densities have been obtained; for the integrated 
partial widths, exact formulae have been derived in the case where the fermions in the 
final state are massless [in Refs ||^, 0], the fermion mass effects have not been taken into 
account and the totally integrated partial widths have not been derived]. The lengthy 
and cumbersome analytical expressions will be found in Ref. [j^. In this section, we will 
simply exhibit the behavior of the x^ and xt branching ratios in the case of decays into bb 
and TV final states, respectively [which are of more immediate interest, being in the range 
probed at LEP2 and the Tevatron]. They are shown in Figs. 4 and 5 for two values of 
tan/3 and the pseudoscalar Higgs boson mass Ma] M2 is fixed to M2 = 140 GeV and we 
assumed universal gaugino masses at Mqut [this leads for /x ~ 450 GeV and tan /3 = 50 to 
the set of gaugino masses: m^o ~ 60 GeV and m^o ~ 171 + ~ 135 GeV: the variation with 
tan/3 is mild]. The trilinear couplings Ai are fixed to 100 GeV for squarks and sleptons. 

In Fig. 4, BR(x2 Xi^^) is shown as a function of (left) and n (right). For large 
values of tan (3, the Abb and hbb couplings are strongly enhanced, and for M4 = 100 GeV, 
the Higgs exchange contributions are largely dominant pushing the branching ratio from 
~ 20% [when only the Z-boson exchange contribution is important] to values close to 
90% [the remaining 10% are in general taken by the branching ratio of the decays into 
r lepton pairs]. For Ma = 500 GeV, the Higgs exchange contribution is suppressed but 
for relatively small values of m^^ or large values of the parameter fi [leading again to a 
small mjj, BR(x2 ~^ Xi^b) is still enhanced. Note that even for moderate tan/?, ~ 5, the 
branching ratio can be large if Ma is small and/or n is large. 
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Figure 4: The branching ratio BR(x2 xV^) ^ function of mi for /j = 450 GeV (left) 
and as a function of fi (right) [assuming in this case, a common scalar mass for squarks 
and sleptons, rh = 380 GeV] for two values of tan /5 = 5, 50 and Ma = 100 GeV (full 
hues) and 500 GeV (dotted lines); M2 is fixed to M2 = 140 GeV. 




Figure 5: The branching ratio BR(xi" — > Xi^'t'^) as a function of m^j for /j, — 450 GeV 
(left) and as a function of /j (right) [assuming, in this case, rfiq = 400 GcV and rhi = 250 
GeV] for two values of tan/i = 5, 50 and Ma = 100 GeV (full hues) and 500 GeV (dotted 
lines); M2 is fixed to M2 = 140 GeV. 
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Fig. 5 displays the branching ratio for the hghtest chargino decay Xi Xi^^^r as a 
function of rrif^ (left) and /i (right). The figure shows the same pattern as for the previous 
decay: for large values of tan (3, a relatively light charged Higgs contribution can strongly 
enhance BR(x]'' —>■ Xi'^'^^r), but even with a heavier boson, the branching ratio can 
exceed the level of 80% because of a lighter fi eigenstate. 



We have developed a fortran code called SDECAY [|rT| which calculates the partial decay 
widths and branching ratios of the chargino, neutralino and gluino decays discussed in this 
noteQ. For the supersymmetric spectrum [including the renormalisation group equations 
for parameter evolution] and for the parameterization of the Higgs sector, it has been 



interfaced with the programs SUSPECT and HDECAY [Q. We have compared our 
results with those of Ref. [^] which have been implemented in the program ISA JET [ p^ . 
For massless fermions, the agreement was perfect]^, giving a great confidence that this 
rather involved calculation is correct. 



4. Conclusions 

In this note, we have analyzed the three-body decay modes of gluinos, squarks and those 
of the lightest charginos and next-to-lightest neutralinos in the MSSM. We have made 
a complete calculation of the decay widths and branching ratios, taking into account 
all possible channels, the mixing in the sfermion sector and the finite masses of the final 
particles, and provided a fortran code for the numerical evaluation of the branching ratios. 

For large values of tan /5, the bottom and tau Yukawa couplings become large, leading 
to smaller masses of the tau slepton and bottom squark compared to their first and second 
generation partners. This leads to enhanced branching ratios of gluinos into bb final states 
and to the possibility that squarks decay into lighter sbottoms through gluino exchange 
can have sizeable branching fractions; this is particularly the case in scenarios where the 
soft SUSY breaking gaugino mass parameters are not unified at the Grand Unification 
scale. The branching ratios of the decays of the lightest chargino into tu final states and 
of the next-to-lightest neutralinos into bb and r"*"r~ pairs are also enhanced in the large 
tan/3 scenario, with an additional increase being due to the stronger Yukawa couplings 
to charged and neutral Higgs bosons, respectively. 

Thus, SUSY events will contain more 6-quarks and r leptons in the final state than 
initially expected. This renders the search for SUSY particles and the measurement of the 
SUSY parameters, where the electron and muon channels where used, less straightforward 
as already discussed in Ref. [|r^ . 6-tagging and the identification of the decays of the tau 
leptons become then a crucial issue in the search and the study of the properties of these 
particles, in particular at hadron colliders such as the Tevatron and LHC 

^The program contains, in addition, the branching ratios for the four-body decay modes [which account 
also for the three-body decays] of the hghtest top squark [|l2| . 

^The comparison was shghtly tricky since the evolution of the couplings and the soft SUSY-breaking 
terms as well as the parameterization of the Higgs sectors are given in different approximations in the 
programs SUSPECT and ISA JET and we needed to use the same input parameters at low energy in both 
programs. We thank Laurent Duflot from ALEPH for his help with this comparison. An independent 
check in the case of the chargino decays, has also been performed by F. Boudjema and V. Lafage n&. 
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